THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200.20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


>JNT  OR  C/TRER  DJIA  WINGS,  SPECIFICATIONS  OR  OTHER  DATA 
; IE  OTHER  THAN  ilN  CONNECTION  WITH  A DEFINITELY  RELATED 
vNT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
-FTY  OBLIGATION  WHATSOEVER;  ANT  THE  FACT  THAT  THE 
FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
ITIONS,  OR  01  HER  DATA  IS  NOT  TO  BE  RlilGARDED  BY 
tSE  AS  IH  ANY  MANNER  UCENSING  THE  HOLDER  OR  ANY  OTHER 
. OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
: £D  INVENTION  THAT  MAY  IN  ANt  WAY  BE  RELATED  THERETO. 


^ V 


J 


0.  N.  R.  EE,SE.\RGH  MEMORANDUM  No.  10 


COMFJTATIONAL  THEORY  OP  LINEAR  PROC-R.'\MmING 


I:  THE  "BOUNDED  VARIABLES"  PROBLEM 


^7 

A.  Chames  and  C.  E.  Lemke 


January  7,  1954 

Graduate  School  of  Industrial  Administration 
Carnegie  Institute  of  Technology 


This  paper  was  written  as  part  of  the  project,  "The  Planning 
and  Control  of  lindustrial  Operations"  under  grant  from  the 
Office  of  Naval  Research.  Tlie  paper  is  released  for  limited 
circulation  in  order  to  facilitate  discussion,  and  Invite 
criticism.  It  is  wholly  tentative  in  character  and  should 
not  be  referred  to  in  publication  without  the  permission  of 
the  authors. 


CQI<tPUTATION;LL  THEORY  OF  LINEAR  PROGRAMMIKa  I; 


THE  «'BGUWDED  VARIABLES"  PROBLEM 
Introduction 

This  paper  is  the  first  of  a series^/  devoted  to  extension  and  sharp- 
ening of  the  present  methods  of  linear  programming  so  as  to  make  practic- 
able the  calculation  of  programs  invol‘';ing  many  more  restrictions  than 
can  now  be  conveniently  handled.  These  extensions  are  developments  on  the 
technique  first  employed  by  Chames  (ref.  [1])  in  reducing  every  linear 
programming  problem  to  one  with  a bounded  set  of  solutions. 

We  here  develop  an  "extended"  simplex  method  for  the  "bounded  variables" 
problem,  e.  g.,  any  linear  programming  problem  in  which  every  variable  en- 
tering Is  constrained  to  lie  between  an  upper  and  a lower  bound.  This  in- 
cludes as  a special  case  those  problems  in  which  not  all  are  so  constrained 
since  we  can  prescribe  arbitrarily  large  (or  small)  upper  (or  lower)  bounds 
for  those  unconstrained.  The  advantage  achieved  v/ith  the  extended  method 
is  reduction  in  size  of  the  computational  tableau  by  suppression  of  the  in- 

zl 

equalities  expressing  the  bounds."' 

Examples  of  such  constraints  occur  repeatedly  in  policy  limitations  on 

items  of  inventory,  in  market  limitations  of  saleable  amounts  of  products,  in 

production  and  delivery  requii’ements  (cf.  ref.  [2]),  etc.  Rather  than  expound 

on  such  examples  here  (which  require  individual  development  for  adequate 

coverage)  we  shall  attempt  to  present,  in  each  paper,  an  example  seemingly 

radically  divorced  from  production  or  industrial  connotations  in  order  to 

gain  one  of  the  most  important  by-products  of  mathematical  formalism,  namely, 

\/  We  wish  to  thank  W.  W.  Cooper  for  suggestions,  stimulation  and  encour- 
agement of  this  series. 

<J  We  iinderstand  that  G.  B.  Dantzig  and  associates  at  Rand  Corporation  had 
developed  a similar  method  (suggested  by  Chames,  Ref.  [1])  in  connection 
with  ’unspecified,  security-restricted  pr-obiems.  Their  rss’ilts  aa  yet 
unavailable . 
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the  recojmition  of  an  analogous  problem  in  a completely  different  field 
from  which  new  insights  n>ay  be  gained  into  the  ori.ginal  problem-  The 
example  presented  in  this  paper  is  that  of  plastic  collapse  of  structures 
(frames). 

We  first  review  the  simplex  method  in  terms  suitable  for  either  the 
’’noimial"  procedure  or  the  ‘'modified"  procedure  or  the  "dual"  method  (ref. 
(33,  M and  [5]). 

The  Simplex  Method. 

The  general  linear  programming  problem  may  be  put  into  the  following 

form,  which  we  shall  refer  to  throughout  as  the  "simplex"  problem: 

n 


1> 


It  will  be  convenient-  both  for  the  review  of  the  simplex  method  and 
for  some  simple  proofs,  to  restate  the  problem  in  vector  notation  as: 


Maximize: 

Zo 

= TI  r*  .n  - 

- -Tt-T 
\f  — 

n 

where : 

= , i - 

and: 

^ ^ ^ “ X^  u* 

[1] 


Maximize: 


where: 


n 

= E 


z-  p.c.  , 

.1=1 


n 


E p.F  , and  Pa  > 0 > j “ n.. 

.1=1  ^ J ^ 


^ .^11  A:i(.4\i.*.i..A»NA9Gua  s/o. 

nucl  matrices.  Throughout  we  shall  refer  to  cj  as  "the  scalar  corresponding 
to  the  vector 

Suppose  that  a set  of  linearly  independent  vectors,  call  them  a^^,  a2» 
...,  are  selected  from  the  vectors  P^,  Pp.  Thus,  Bj^  is  a basis  of  Vj^ 

and  we  may  vrrite  tiniquely: 

m 

Pq  = 2 P^  • 

i=l 


[2] 


If  the  m values  satisfy  > 0;  i = !)  , thr-;'-  th'-..  "oi-m  a soJ.ution 

to  the  system  [1].  iifhen  the  also  satisfy  p^  > 0*  i = 1,  ....  m,  the  solutiuu 
is  referred  to  as  a basic  solution.  It  has  been  shown  (ref.  [1])  that  the 
problem  ca)s  be  so  modified  that  when  one  has  a solution  to  [1]  of  the  fora 
[2]  the  soj.ution  is  a‘'v;ayG  a basic  solution.  Hence  we  may  make  the  initial 
assumption  that  the  original  set  of  vectors  P^,  P^,  Pjj  has  this  property, 

wiiich  may  be  characterized  by  the  statement  that  the  vector  ?q  is  linearly 
independent  of  any  m~l  vectors  chosen  from  among  Pq^,  ...,  P^. 

With  this  assumptj.on,  the  simplex  method  may  be  described  as  one  which 
proceeds  from  basic  solution  to  basic  solution  until  a maximal  basic  solution 
is  found,  i.  e.,  one  which  yields  the  maximum  value  for  Zq.  We  shall  review 
the  procedure  for  going  from  one  basic  solution  to  another  by  the  simplex 
method.  Thus,  suppose  one  has  located  a basis  Ej^  such  that,  expressing  Pq 
as  in  [2]  one  has  pj[  > Oj  i = 1,  ...,  m.  Let  be  the  scalar  correspon- 
ding to  the  vector 

We  introduce  the  ’vinique  vectors  a^j  i = l,  ...,m  satisi.ying: 

[0  if  i j 


[3] 


^i 


*a^  **  ^ii  " ^ i-or  ^ , J — 1,  • • • , m. 

[1  if  i = J 


We  shall  call  the  vectors  a^  the  "dual  vectors"  for  the  set  a^,!  If  A is  the 
matrix  whose  ith  column  is  aj_  then  the  jth  row  of  A ^ is  a^.  Having  the  dual 
vectors  and  denoting  bhe  inner  product  of  vectors  x and  y by  x’y  any  vector 


X in  can  be  expressed  by: 


C4] 


rri 

y j 


[6] 


^ i- 

X « S (x*a  )a^ , 
i=l 

In  particular,  we  have:; 

i»*l 

El 


P.  =«  E (p4’a'^)a^  ; j = 1,  ...,  n. 

J ^ -I 
x“X 


Having  the  compated  data  Pq’^  and  Pj'a'' 
ced^are  is  to  corapate  the  rniantities; 

ni 


. the  first  step  in  the  pro- 


[7] 


:(Pj)  = (Pj'a^)ci.^  “ j - 1,  n. 


Z(Pj)  is  formed  b;y  replacing  each  a.j  by  its  corresponding  scalar  Cj.^ 
in  the  expression  [6]  for  Pj  and  then  subtracting  the  scalar  corresponding 


to  P., 


When  the  n quantities  [7]  are  all  non~negative,  one  may  demonstrate 
that  the  pronle™  is  finished:  the  basic  solution  is  a maximal  one.  When  at 
least  one  of  the  quantities  [7]  is  negative,  a new  basic  solution  may  be 
secured  which  increases  the  functional  Z^,  This  is  accomplished  by  forming 
a new  basis  from  3^^^  by  replacing  one  of  the  af  by  another  vector.  The 
replacing  vector  my  be  any  Pj^  for  vAiich  Z(r^)  < 0,  To  decide  upon  the  vec~ 
tor  to  be  replaced  one  finds  the  value  of  i,  say  r,  such  that: 


lS]  © « CPo*an/(Pk*a^')  - min;(Po»a^)/(Pj^’a^)  for  P^*a^  > 0. 

i 

The  following  facts  may  then  be  demonstrated  (ref.  [6l)j 

i;  The  value  r which  yields  9 is  unique, 

ii)  The  set  3«’  consisting  of  Pj^,  ®^r+l>  ****  ^ ^ basis 

and  again  yields  a basic  solution, 

iii)  The  new  basic  solution  yields  a larger  valiie  of  the  '’u'  ctional  Zq, 

Having  decided  upon  the  values  of  k and  r,  when  the  vectors  Pq,  ... 
Pjj  are  expressed  in  terms  of  the  new  basis  the  transition  to  the  new 

basic  solution  may  be  considered  co.mplete.  (There  is  a simple  algorithm  for 
this  transition  which  vd.ll  be  presented  later). 

One  may  then  show  that  in  a finite  number  of  steps  one  will  obtain  a 


maximal  basic  solution,  namely  one  for  which  all  of  the  quantities  Z(P4) 
are  non -negative. 
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It  is  usual,  both  for  computational  and  explanatory  purposes,  to  con- 
struct a oomputational  tableau  which  exhibits  the  computations  pertaining 
a given  stage  and  at  the  ssLine  time  clarifies,  for  the  computer,  the  mechanics 
of  proceeding  to  the  next  stage. 

Computational  Tableau 


Corresponding 

sc» 

1 

alars 

^ 

Cl  ... 

. . . Cjj 

i 

Bm 

Pc 

Pi  ... 

Pk 

• • • Pn 

cri 

• 

"1 

ai 

• 

P*  * 

O 

• 

Pi 'a-  . . . 

• 

• 

P *al 

...  a 

• 

• 

• 

• 

a_ 

Po*a® 

• 

« 

Pl*a3  ... 

• 

• 

XU 

• 

. . . Pn'as 

• 

• 

• 

• 

• 

« 

• 

Pl^a®  ... 

• 

Va® 

• 

. . . Pn*a® 

n 

"o 

Z(Pl)  . . . 

z(pp 

. . . Z(Pn) 

We  shall  refer  to  such  a tableau  as  an  ’’mxn  tableau"  although  there  are 
a few  additional  rows  and  columns,  and  shall  likewise  refer  to  the  problem 
as  an  "mxn  prob] em" . 

The  Bounded  Variables  Problem 

First  an  illustration.  The  following  problem  arises  in  the  thecr.'  of 
plastic  collapse  for  structures  (refs.  [?]>  [8]). 

Consider  a structure  composed  of  elastic,  perfectly  plastic  beams 
joined  together  rigidly.  V7e  are  interested  in  plastic  collapse  of  this  struc- 
ture under  the  action  of  concentrated  applied  forces.  We  assume  that  each 
beam  is  homogeneous  so  that  the  same  yield  conditions _apply  to  any  cross- 
section  of  the  same  beam. 

Consider  a sot  of  numbers  p]^,  ...,  py^  vdiich  represent  the  magnitudes 
of  extemal  forces  which  are  to  be  applied  to  the  structure.  The  structure 
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will  be  in  equilibiriiim  when  these  forces  are  balanced  by  the  reactions  at 
the  supports.  These  reactions,  together  with  the  applied  forces,  give  rise 
to  a system  of  internal  fcroes  which  determine  the  bending  moment  at  every 
point  of  the  system  (the  effects  of  shear  and  asdal  forces  can  be  neglected 
here) . 


For  a given  Cinss-section  there  will  be  a maximum  and  a minimum  bending 
moment  which  the  cross-section  can  carry  due  to  the  assuiaption  of  perfect 
plasticity.  Thus,  at  each  cross-section  the  bending  moment  M vdll  be  subject 
to  an  inequality  < M < H°.  is  the  '’fully  plastic  moment”  or  the  mo- 

ment for  which  plastic  flow  v/ould  first  begin.  Now,  the  bending  moment  dis- 
tribution will  be  linear  along  each  member.  Hence,  we  need  only  be  concerned 
with  those  points  of  the  structure  vdiere  plastic  yielding  will  first  occur. 
These  points  are  finite  in  number  and  depend  only  on  the  structure  and  on 
the  points  of  application  of  the  forces.  They  are  points  including  those 
for  which  the  bending  moment  has  a turning  point,  as  well  as  the  ends  of  the 
ueams.  Ox*  whex*e  a support  3.3  present,  or  poants  at  whach  one  of  the 

forces  is  applied.  These  points  are  called  "test  stations”.  Thus,  at  each 
test  station  the  bending  moment  Mj  is  subject-  to  the  restrictions: 


C9] 


o x> 

-M.  < hj  Pij 


i j ^ j ^ j } j — 1,  ..  .,n. 

Now,  considering  the  forces  p]^,  . . . , p^^^  as  a base  load,  v.*e  apply  instead 
the  loads  pp]^,  ...,  ppi^.  Then,  besides  the  restrictions  [9],  the  bending 
moments  will  bo  subject  to  the  equilibrium  conditions: 
n 

[xOj  E bj  jMj  — pp.  j X = 1,  ...,  m, 

j=l  “ 

vdiere  the  b^^  depend  on  the  structure  and  on  the  points  of  appli  cat X on  of 
the  forces. 

uC  scale  the  largest  value  P of  p permxssxbls  xn  order  tnau  the  structure 


remain  in  equilibrium  under  tne  action  of  the  loads  pp^ 


The  prob- 


We  make  c.he  following  tran;:f'onr.2.tions  and  nctalional  changes: 

n 

'^1  = P ; 3 - . = ; aiid  E a^j. 


Let  Xj  - + 1 


lem  then  becomes: 


where: 


and: 


^+1 
n 


0 ^ X j ^5  J 

This  is  an  example  of  the  follov/ing  problem  which  we  shall  state  in 
vector  form  and  which  we  shall  call  the  ”bounded  variables"  problem: 


Maximise: 

where: 


2 Pj<^i  > 
1=1 


n 


Z ^ ‘ J ~ • » m f TX  a 


W'j  sijiall  set  up  and  treat  this  problem  as  a simplex  problem.  In  this 
manner  x-rs  shall  have  an  (m+n)x(2n)  problem.  Such  a problem  would  ordinarily 
employ  (m+n)x(2n)  tableaux.  We  shall  show,  however,  that  it  will  be  suffi- 
cient to  employ  roxn.  tableaux;  that  is,  that  the  problem  can  be  handled  in  the 
space  Vm  of  the  vectors  Pj, 

If  we  introduce  the  new  variables  x^  defined  by: 
tllj  Pj  Xj  = bj  > j=l,  ,,,,  n, 

wa  may  phrase  the  problem  as  follows  in  simplex  form: 


Maximize: 


wnere: 


O — 

n 


Po  = 2 PiP.  » 
j=l  ~ 


’j  - pj 


fU  * X. 


•j 


and 


Pj>  Xj  > 0 ; J =•  1, 


n. 


We  wish  to  define  vectors  in  These  will  be  designated  by  putting 

_ _ Tx? 

bars  above  tho  letters,  as  z.  At  times  we  shall  use  the  notation  z = 

where  x is  a vector  in  and  y is  a vector  in  V.,, 


e " 


Vie  introduce  the  following  vectors  in  V_: 

n 


D “• 


/x-.\ 

“1 


bn 


l^\ 

'■•i 

n I 


\°l 


-(;3th  place)  j j = 1,  n. 


And  the  following  vectors  in 


rpi 


= 


J 

[Qj 


■ 


■°1 . 

L«jj 


; j = 


Tiie  (in^^n)x(2n)  prcblssi  may  then  be  phra‘’e4  as: 

n 

Maximize:  2q  - E pjCj  , 

r 12  3 an 


n. 


vrtiere : 


E PjPj  + E X 


J=1 


j=i 


^ ^ > 3 ~ l>t»»>n» 


Suppose  that  a basic  solution  for  this  pi’oblem  has  been  fovnd.  This 
implies  that  a basis  B^+jj  of  has  been  chosen  from  among  the  vectors 

P^,  Pnj  Qi  > •-  = , Qn«  Recall  that  the  simplex  method  req’iires  that  each 

of  these  2n  vectors,  together  with  Po,  be  expressed  in  terms  of  the  basis 

We  show  that  the  coefficients  involved  can  be  obtained  by  merely  ex*- 

m+n 

pressing  the  vectors  P]^,  Pj^  in  terms  of  a certain  basis  of  Vm. 

To  do  this  let  us  first  note  some  properties  of  an  arbitrary  basis 
E-+V,  chosen  as  above- 

i)  The  number  s of  vectors  Pj  in  Bjjj+jj  satisfies  m < s < n. 

ii)  For  each  j,  either  Pj  or  Qj  or  both  are  in  ®m-«-n* 

iii)  By  i and  ii  there  are  precisely  m values  of  j such  that  both  Pj 
and  are  in  B^j+jj*  Furthermore,  the  set  B^  of  the  m vectors  Pj  such 
that  both  Pj  and  Qj  are  in  is  a basis  of 

In  particular,  suppose  that  B^^+jj  yields  a basic  solution  for  the  prob- 
lem [12].  For  notational  ease  we  suppose  that  I^,  Pg  are  in  Bja+yj  and 
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( ; 


that  and  are  both  in  B^+n  -^or  i = 1 , . . . , m.  Asa  result  we  have  that 
Qj.  complete  the  basis  3^^+-  and  that  ^ > Pjj^  form  a basis 

of  Vjj. 

Let  a^.,  a“  be  the  vectors  dual  to  the  basis  Thus,  we  have; 

[133  = Oij  ; i>  J = 1.  m. 

Suppose  that  we  have  computed  the  quantities  Po*a-^  and  Pj*a^  satisfy- 
ing; 


m 


riM 

L j 


L15; 


= E (Po'al)Pi 

i=l 


Pj  « 2 (Pj*a^)Pi  I 1 = 1,  ...,  n. 


By  means  of  these  quantities  we  may  write  down  the  expressions  for 
"P^,  "pj  for  s < .j  < n and  Qj  for-  m < J < s in  tenns  of  3,n+n!  namely,  for 
those  vectors  not  already  in  the  basis  Thus,  we  have: 

m 


[15]  Pj  “ Qj  + 

[16]  Q. 


w - 
o 


1P4 

J 

IDJ 

rv>  ^ 


[17] 


P =1 

o 


» Qj  + E (p4»a^)rPi  - Q^]  ; s < j < n, 

J i=l 

= Pj  - 2 (Fy%^)[F^  : m < j < s. 


i=l 


ID  J 

'0  ~!  m „ _ n _ 


+ I I ” 2 ( 

Lb  J i=l 


m 4_«_  ® “ • 

2 (Po’a^)[Pi-Qi]  + 2 b^CL  - E b^Q^  + E b rp^  - E P^ »a^(Pi*^i) ] 
i=l  ^ i=l  i=s+l  ^ j+m+1  ^ ^ i=l 

- E(5'iPi+  E(bi-0^)C^  ' 

i=l  i=l  J"mtl  J=s+1 


•>  ^ • 

where  we  have  set  0^  * Po*a^  E b.(P4*a^)  for  i = 1,  . . . , m. 

J=m+1  '' 

We  recall  now  that  the  first  step  in  applying  the  simplex  method  is 
to  calculate  the  quantities  2(Pj)  and  Z(Qj),  obtained  (see  eqn.  [?])  from 
the  vector  expressions  [15]  and  [16]  by  replacing  each  of  the  baeds  vectors 
by  its  corresponding  scalar. 


Thus,  in  view  of  [15  3 Olid 


[16],  these  quantities  become: 


[18] 

Z(pj)  « 0 + 

[19] 

z(Qj-;  . c,- 

m 


m 


I’a  )c^  - » Z(P.)  : s < i <r  n. 

I j.  j j ' ~ - 

Z(Qj)  = c^-  S (Pj*a^)ci  - 0 = -Z(Pj)  ; m<  j < 3. 

The  values  in  [IS]  and  [19]  will  determine  whether  or  not  the  basic 
solution  corresponding  to  the  basis  is  maximal.  If  it  is  not,  one  of 

these  quantities  will  be  negative  and  we  may  proceed  to  a new  basic  solu- 
tion. In  this  case,  having  decided  which  vector  will  replace  a vector  in 
W0  ri6sd.  1jo  ^*ind.  ths  vsctoi*  v;'';\ch  is  to  be  I'oplscod*  Foi*  this  hs 
need  the  quantity  Q (see  eqn.  [6]).  The  ratios  competing  for  Q are  found 
by  means  of  equations  [15j»  [16],  and  [17]»  We  need  to  distinguish  two 
cases, 

i)  If  Qj^  for  m < k < s is  to  replace  some  vector  in 
© » Min.  f[Min.  0i/(-Pv»a^)  for  Pi-’a^CO],  [Min. (b, )/Fvf for  Pic’ai>0- 

X 1 J 

ii)  If  Pjj  for  s < k < n is  to  replace  some  vector  in  EU+nS 
9 = Mn.  |"[Min.0i/Pj^’a^  for  Pj^'a^>0],  [Min,  (bj^-04  )/(-Pj^’a^)  for  I^’a^^],  b^  ^ 

t.  i i J 


Having  determined  the  vector  to  be  replaced,  one  can  complete  the  re- 
placement and  proceed  to  the  new  basic  solution. 

We  have  thus  shoim  that  essentially  all  of  the  computations  in  one  stage 
of  the  (m+n)x(2n)  problem  may  be  handled  by  means  of  an  mxn  tableau,  Tr  com- 
plete the  discussion  of  the  bounded  variables  prt^blem,  we  shall  exh’.bit  the 
computational  tableau  associated  with  one  stage  of  the  problem  and  list  the 
rules  of  procedure  tc  be  followed  in  passing  to  the  next  stage. 
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Explanations  and  Ptules  of  Frocedure- 

1,  The  basis  is  vinioalized  in  the  tableau  as  follows: 

Pj  is  in  if  the  last  entry  in  the  column  headed  by  Pj  is  a +. 

Qj  is  in  if  the  last  entry  in  the  column  headed  by  Pj  is  a 

Thus,  both  + and  - at  the  bottom  of  the  column  headed  by  Pj  indicates 
that  both  Pj  and  Qj  are  in 

For  the  procedure,  one  first  cheeks  the  n values  Z(Pj)  in  the  next  to  the 
last  row,  A maximum  has  been  reached  if  both: 

Z(Pj)  > 0 in  all  of  those  columns  having  + in  the  last  entry; 

Z(Pj)  < 0 in  all  of  those  columns  having  - in  the  last  entry. 

[This  would  mean  that  all  of  the  quantities  Z(Pj)  and  Z(Qj)  were  non-negative,] 
If  some  Z(Pj)  violates  the  above,  we  proceed  to  a new  basic  solution. 

Having  decided  upon  the  replacing  vector,  v'c  need  the  value  0 given  in  either 
i or  ii  on  the  preceding  page.  It  is  for  this  reason  that  the  columns  a and 
b under  Po  and  the  first  row  were  added  to  the  tableau. 

Having  decided  upon  the  replacing  vector  and  the  vector  to  be  rsplaced 
may  pass  on  to  the  new  basic  sclut.ion  and  the  new  tableau.  The  following 
cases  may  arise: 
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i*  The  vectji'  replaces  the  vector  Pj^  in  Sjjj+jj  Pcx"  ijoiue  k with  in  < k ^ s. 
This  induces  no  change  in  the  basis  The  only  changes  in  the  tableau  in 

this  case  are  the  foliovd-ng: 

The  + under  hhe  column  headed  by  Pjj.  !•?  changed  to 
The  quantity  0^*  « 0j_  + bjfCPj^’a^) : i = 1,  ..  .,m  replaces  0^  in  the 
a and  b columns  under  Pq.  (see  the  expression  for  0^  in  oq'*.„  I7). 

_o.  The  quantity  bi^Z(Pj^)  ir-  \,o  Z^.  (Z^,,.  of  course,  need  not  be 

carried  along.  The  maximum  value  of  Zq  may  be  obtained  at  the  end  of  the 
problem  from  eqn.  17 •) 

ii.  The  vector  Pi,  replaces  the  vector  Q<„-  in  for  soiae  k with  s < k < n. 

■»  1^  XV  “ xw  *IA.  • *X  ““ 

•'his  induces  no  change  in  The  only  changes  in  the  tableau  in  this  case 

are  the  following: 

_a.  The  - under  the  column  headed  by  Pv  is  changed  to  +. 
b«  The  quantity  0^*  = 0i  - b»^(Pj^*a^)j  i = 1,  ,..,m  replaces  in  the 
a and  b columns  under  p^, 

jc.  The  quantity  bi^ZCPk-)  is  subtracted  from  Zq. 
ill.  If  either  Pj.  or  Q^.  with  1 < r < m is  replaced  either  by  some  P^-  with 
s < k < n or  some  Qj^  with  m < k < s then  replaces  Pj.  in  the  basis  Bm  to 

form  a new  basis  Bjjj*  anu.,  as  in  the  usual  siiiipleji.  pr-oCeuui-c,  the  I'lcW  tableau 
must  be  based  upon  the  basis  Bj- ' . The  simplex  algorithm  for  changing  the 
tableau  elements,  which  is  included  here  for  t’ne  sake  of  coinple  ceness,  is 
as  follows.  The  element  is  replaced  by  the  element: 

(Pj’a^)*  = Pj’a^  - (P[^*a^/Pif ’a")Pj'a^  for  j = 1,  n but  j / k, 

and  for  i=l,  m but  i / r, 


(Pj'a^)*  - Pj»aV>^'a“ 

^ [0  if  i 0 k 

(P^»ai)»  » 


;;  j = 1, 


r 1 • ^ j 


The  n quantities  Z(Pj)  are  replaced  by: 

Z(Pj)»  = Z(P^)  - [Pj*a7Pk’a^3-'^<Pic)  5 J = 1,  n. 
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Und<=ir  pQ  the  nol’amne  headed  a arid  b are  changed  by  replacing  by: 

* = 0^  - ; Tor  1 = 1,  m,  but  i f r,  and  €»  replaces 

0J.  in  the  rth  place  of  the  a eolaiuDj  - S replaces  bj,  - 0^.  in  the  rth 
place  of  the  b column. 

Finally,  replaces  Pj.  in  the  column;  replaces  c,..  in  the  corres* 
ponding  scalars  column;  and  both  + and  - appear  in  the  last  entrj’-  of  the 
column,  whereas  the  P-  mliirnn  ha«  a + if  Qj.  has  been  replaced  in  or  a 

- if  Pj.  has  been  replaced. 

This  completes  the  discussion  of  the  bounded  variables  problem. 

As  a final  remark,  it  should  be  noted  that  when  a ctiange  in  the  tableau 
(case  iii  above)  is  necessary,  the  modified  simj)lex  method,  (ref.  [4])  may 
be  employed,  whereby  only  the  vectors  (a^)  * dual  to  the  new  basis  Bm*  are 
computed  from  the  old  vectors  a^,  and  then  the  quantities  (Pj'a^)*  are  ob- 
tained, if  needed,  by  taking  the  scalar  product  of  Pj  and  (a^)’.  The  al- 
go H thm  for  obtaining  the  new  dual  vectors  is: 

(a^)*  = 

(a^)*  = -(Pj^’a^)  (a^)  * for  i = 1,  m but  i / r. 
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